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Ocean Floor Topography 





Oceanographer Marie Tharp is shown in this 
photo from the early 1950’s at her drawing board de- 
tailing one of her famous geophysical charts of the sea 
floor. She has been associated with Lamont-Doherty 
Geological Observatory, Columbia University, since its 
founding in 1949 until retiring in 1983. During this 
time she collaborated with the late Dr. Bruce C. 
Heezen in the study of sea floor topography; this work 
was funded by the Office of Naval Research for many 


years. In 1977 their studies resulted in the famous 
World Ocean Floor Panorama. This pictorial represen- 
tation of the sea floor is based on more than 5,000,000 
miles of soundings taken by oceanographic institutions 
and hydrographic bureaus of many countries and the 
bathymetric studies by Dr. Heezen and Marie Tharp. 
The chart remains today a standard for its accuracy in 
depicting the features of the ocean floor. 
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Figure 2: 


Sierpinski gasket is a regular fractal made up of triangles. 
One starts with a single triangle. In the next step, three 
triangles are put together to form a new triangle with a 
triangle missing in the middle. In the next generations, 
again three of the new triangles are put together to form a 
larger triangle with a triangle missing in the middle. Clearly 
this process can be continued ad infinitum to form a very 
large pattern. Since the repeating unit is always of the form 
of three triangles with a missing triangle in the middle, this 
object is self-similar at any length scale. The fractal 
dimension of a Sierpinski gasket is D=log3/log2=1.58. . . . 
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In contrast, fractal shapes and patterns found in nature are 
usually random fractals. The reason is that they consist of 
random shapes or patterns that appear stochastically at any 
iength scale. Because of this randomness their self-similarity 
is only statistical. Another property of natural fractals is that 
their self-similarity extends only over a finite range, for ex- 
ample from atomic sizes to the maximum length in the sys- 
tem. Examples of random fractals include such objects as the 
coastlines and mountains, mentioned earlier as well as many 
disordered materials like polymers, colloids and gels. In fact, 
much of the current activity in the application of fractal con- 
cepts to physical systems has been in the area of random 
materials.* 

The fractal dimension of an object or a pattern serves 
to characterize it and often indicates the manner in which it 
has been formed. There are several ways to determine the 
fractal dimension, but one of the most widely used methods 
is to relate the size of the object to its length. Let us consider 
the Sierpinski gasket shown in Figure 2. In each generation 
of the Sierpinski gasket, the length of the new triangle is 2 
times the length of the previous one, and its mass is 3 times 
the previous mass. In contrast, if we construct a non-fractal 
object, like a regular triangle, then every time its length is 
doubled, its mass is quadrupled. If the length of a regular 
triangle is increased by a factor of n, then its mass would 
be n?M(r), where r is the length of the original triangle and 
M(r) is its mass. The number 2 is the dimension of the ob- 
ject; a trivial result for a regular two dimensional triangle. 


In general, the fractal dimension is the exponent or power 
to which n is raised. The fractal dimension D is defined by 
the relation, M(nr)=n?M(r), where M(nr) is the mass after 
the length has been scaled by a factor of n. For the Sierpin- 
ski gasket M(2r)=3M(r), because every time r is doubled, 
the mass is tripled. Therefore, the fractal dimension D of a 
Sierpinski gasket can be determined from the relation, 
3=22, which implies that D=log3/log2=1.58.... This 
method can be used to determine the fractal dimension of 
any object or pattern. In the case of a random fractal, one 
can simply draw circles of increasing radii about a point on 
the object and determine the mass M(r) within a circle of 
radius r and compare it with the mass within a circle of radius 
nr. The fractal dimension D can be calculated from the rela- 
tion M(nr)=n?M(r), as before. 

Time can also be fractal. Let us assume that we use a 
certain fixed time interval to make a standard measurement 
of time. This could be, for example, beeps recorded on a mag- 
netic tape every one second. We can tell time by this device 
no matter how fast or slow the tape is played. If the tape is 
played at twice the speed, the time between beeps is half a 
second; if it is played at half the speed, the time interval is 
two seconds, and so on. Now, imagine a process consisting 
of events or beeps that is fractal in time. In such a process 
there is no characteristic time between beeps. If we make 
a recording of the events in this process and we play the tape 
at any speed, we find beeps at all possible time intervals: 
from the shortest time that we can measure between two beeps 
to as long as we listen to the tape. There is no fixed or charac- 
teristic time interval that we could use as a standard to meas- 
ure the time between beeps. We cannot use the average time 
between beeps either, because it depends on the length of time 
that we listen to the tape. This is the distinguishing feature 
of a distribution of times that is fractal. The reason is that 
for a random distribution of time intervals, the average time 
is always constant. No matter how long we listen to the tape, 
the mean time between the beeps is the same. Processes with 
a fractal time distribution occur commonly in nature; the best 
example is 1/f noise. 

In the following, I will briefly discuss some of my re- 
cent work on the development and application of fractal con- 
cepts in physical systems which have been supported by the 
Office of Naval Research. Despite the limited scope of this 
review, I think it should give an indication of the richness 
and complexity of the phenomena that can be understood 
using the mathematical tools of fractal geometry. 
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Percolation and Gelation 


Many materials of technological and practical importance 
have an inhomogeneous, disordered structure somewhere be- 
tween ordered crystals and random gases. Examples include 
the gelation of polymeric materials* and the growth of 
clusters on surfaces via atomic chemisorption, composite 
mixtures of conducting and dielectric materials, and dilute 
magnets. The formation of this type of system can be 
described by a growth process in which individual elements 
randomly join together to form an interconnected, percolat- 
ing network.** The most widely used model for such sys- 
tems is the percolation model, in which sites (or bonds) on 
a lattice are occupied randomly with a probability p, indepen- 
dent of the occupation of the neighboring sites. Nearest- 
neighbor sites are then considered to be connected, result- 
ing in clusters of varying sizes and geometries. If the proba- 
bility p is large enough, p>p,. one cluster will percolate 
through the entire lattice. This “infinite” cluster corresponds 
to a gel in gelation, or to a conduction path through metallic 
atoms in a mixture of conducting and dielectric materials. 
Its formation at the percolation threshld p, signals a geomet- 
ric phase transition, and it is well established that in the re- 
gion around p,. random percolation exhibits scale-invariance 
and fractal behavior. In particular, the percolating cluster at 
the percolation threshold is a fractal object, with holes and 
irregularities at all length scales. 

In general, however, a physical system will have some 
type of interaction between its individual elements.’ In sur- 
face deposition, for example, the presence of neighbors can 
either enhance or deter absorption at a site. Such interactions 
will give rise to correlations which can affect the properties 
of the system in important ways. We have introduced and in- 
vestigated an interactive percolation model which is a sim- 
ple generalization of random percolation.** Instead of using 
a single probability p to decide if a site should be occupied, 
two probabilities are used, pp and p,, according to whether 
the site has no neighbors or at least one neighbor, respec- 
tively. The value of the percolation threshold p, and other 
properties of the system depend on the parameter r=p,/po. 
For r=1, the system corresponds to random percolation. 
When r>1, the site occupation will be “ferromagnetic” in na- 
ture, since sites that have occupied neighbors will be pre- 
ferred over those that do not. Cluster nucleation will be small, 
resulting in widely separated clusters. Eventually, these com- 
pact clusters, or “blobs”, link together and form large, fractal 
clusters, one of which will percolate through the system. This 
can be observed in Figure 3, which shows the distribution 
of clusters at the percolation threshold for several values of 
r. When r<1, on the other hand, the site occupation will be 
“antiferromagnetic” in nature, since the presence of occupied 
neighbors will be inhibitory. Initially, therefore, there will 
be many single-particle blobs, which will only link together 
with difficulty. They will tend to lie on one of two sublat- 
tices, forming next-nearest-neighbor-connected domains sepa- 
rated by antiphase boundaries. 
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At first glance, the percolative structures shown in Fig- 
ure 3 appear to be very different from each other. Consider, 
however, when the value of r is large, and instead of the origi- 
nal particles, the blobs are taken to be the basic elements 
which join together to form percolating clusters. Then, these 
pictures are very similar to random percolation (on a smaller 
lattice). This suggests that as long as all quantities are scaled 
by the linear size of these blobs, the properties of the system 
will not change. We have made extensive investigations of 
the properties of the interacting percolation model. We find 
that our data for all finite r are consistent with the idea that 
the fractual dimension and other properties of the interac- 
tive percolation model are the same as random percolation. 
Thus, the interactions introduce correlations that are signifi- 
cant only on length scales smaller than the blob size. 


Fractal Distribution 
of Clusters 


The formation of a distribution of various size clusters 
is the characteristic feature of numerous processes in physics, 
chemistry, biology, and engineering. Despite the fact that the 
formation of clusters, from atomic clusters to a distribution 
of galaxies, had been a subject of intensive investigations for 
many decades, no unified description of these processes ex- 
isted. Thus, the recent development that the distribution of 
clusters is scale-invariant °* and can be described by frac- 
tal concepts has opened the way for important new develop- 
ments which are leading to a unified description of many 
cluster-growth processes.?:>8 

In general, cluster formation can be divided into two 
classes depending on whether the clusters have a ramified, 
fractal structure or are smooth and spherical like liquid 
droplets. Fractal clusters are generally found in aggregation 
processes *.°'° involving small particles, which include col- 
loids, aerosols, polymers, and many biological processes. In 
contrast, droplet growth '*"'* is observed when the coalesc- 
ing clusters are in the liquid state, such as the formation of 
vapor-deposited thin films, breath figures, clouds, rain, fog, 
foam, and froth. Figure 4 shows an example of a model of 
cluster growth by diffusion limited aggregation ? which leads 
to the formation of a distribution of fractal clusters. In this 
process the system initially consists of a fixed number of - 
monomers which diffuse and interact to form larger clusters. 
Experiments and simulations have shown that large aggregates 
have a fractal geometry with a dimensionality D~1.4 in two 
dimensions and D~-1.8 in three dimensions. 








Figure 3: 


Sample cluster distribution at the percolation threshold in 
the interactive percolation model for various values of r. The 
percolating cluster is black; all others are gray. For r=1, the 
system corresponds to random percolation. When r>1, 
cluster nucleation will be small, resulting in widely 
separated clusters. When r<1, the presence of occupied 
neighbors will be inhibitory. Initially, therefore, there will be 
many single-particle blobs, which will only link together with 
difficulty. They will tend to lie on one of two sublattices, 
forming next-nearest-neighbor-connected domains 
separated by antiphase boundaries. [From ref. 3] 
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(a) r = 0.001 


(f) r = 10,000 


Figure 4: 


Evolution of the cluster size distribution in a diffusion-limited 
aggregation process. The system is initially composed of a 
fixed number of monomers. The monomers and clusters 
diffuse and join together to form larger clusters. The fractal 
dimension of the large clusters in two dimensions is 1.4. 
[From ref. 11] 















































There are many processes in which spherical droplets 
grow by absorption or condensation of smaller droplets or 
atoms.'*"'® An example of the formation of a distribution of 
regular shape droplets is shown in Figure 5, which is an SEM 
micrograph of tin vapor deposited on a sapphire substrate. 
In this type of experiment ' the substrate is held at a tem- 
perature slightly higher than the melting temperature of the 
deposited material. Since the deposited material is in a liq- 
uid state, it forms spherical (or hemispherical) droplets on 
the surface which grow and coalesce to form a distribution 
of droplets. The smooth shape of the droplets is preserved 
throughout the growth process, because when two or more 
droplets touch, they coalesce to form a new droplet with mass 
conservation. An example of distribution of droplets obtained 
in computer simulation of a simple model of droplet growth 
is also shown in Figure 5. The agreement between the results 
of the simulations (right) and the experiments (left) indicates 
that our model captures the essential features of droplet 
growth in thin film experiments. 
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SEM micrograph of tin vapor deposited on a sapphire 
substrate is shown on the right and the distribution of 
droplets obtained by computer simulation of a simple 
model of droplet growth is shown on the left. There are 
many qualitative similarities between these pictures and we 
have also shown quantitatively that the model captures 
many of the essential features of droplet growth in thin film 
experiments. [From ref. 16] 


Figure 6: 


Self-similarity of cluster size distribution in aggregation and 
droplet growth processes. In each graph the distributions 
for various times are found to collapse to a single curve 
when they are scaled according to the dynamic scaling 
theory. This type of self-similarity has been verified in a 
wide variety of theoretical and experimental systems, 
including colloids, aggregates and droplets. [From ref. 9 
and 12] 
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The significant result of these studies is that the distri- 
bution of both fractal clusters in such processes as colloidal 
aggregation and smooth droplets in the growth of thin films 
or breath figures are self-similar and can be described by frac- 
tal concepts.'*"* In particular, using the dynamic scaling the- 
ory,'*'® we have shown that the cluster size distribution is 
self-similar in both time and space. This implies that the clus- 
ter size distributions obtained at different times look the same 
if time and length are properly scaled. For example, Figure 
6 demonstrates the self-similarity of cluster size distribution 
in both aggregation and droplet growth processes. In each 
case the distributions for various times are shown to collapse 
to a single “universal” curve when they are scaled according 
to the dynamic scaling theory. This type of self-similarity has 
been shown to exist in a wide variety of systems, including 
colloids, aggregates and droplets. 


Rough Surfaces: 
Self-affine Fractals 


The formation of complex structures and patterns is com- 
mon in a wide variety of physical, chemical and biological 
processes, including crystal growth, solidification, 
aggregation, vapor deposition, and combustion. The most 
characteristic feature of these patterns is the existence of an 
evolving interface separating the structure from the outside. 
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“Diffusion-Limited Aggregation” 





Ds ws? 


In[s?hg(t)] 


—- nN We Aaoan @ Oo 








°o 


-2 wee 
In (s/S(t)) 


“Droplet Growth by Homogeneous Nucleation” 





T T T 





d=2,D=3 
A= 512 x 512 
Mmax = 1134878 








L lL 
-5 -2 
In (s/S(t)) 








Thus, one way to develop a better understanding of growth 
phenomena would be to study the dynamics of the inter- 
face.'*-*! Surface properties play an important role in scien- 
tific and industrial applications as well. Therefore, charac- 
terization of surface roughness is a problem of considerable 
practical importance. 

We have studied the evolution of the interface in a num- 
ber of models which simulate growth by vapor deposition or 
sedimentation.'**! In this class of models we let particles 
simply rain down and stick to a surface. Newly arriving par- 
ticles stick and become part of the aggregate when they reach 
a particle that is already part of the deposit. Two examples 
of the type of surfaces found in these simulations are shown 
in Figure 7. In the process on the right we have incorporated 
the effects of local diffusion by allowing a newly arriving par- 
ticle to diffuse on the surface until it comes to rest at a local 
surface minima. This effect which tends to smooth the sur- 








Figure 7: 


Two types of deposits and interfaces that are found in the 
simulations of models of thin film growth by vapor 
deposition or sedimentation are shown here. In this class of 
models particles simply rain down and stick to a surface. 
Newly arriving particles stick and become part of the clus- 
ter when they reach a particle that is already part of the 
deposit. In the process on the right particles fall to a 
position of local minima which tends to smooth the surface 
and form a compact deposit. In ballistic deposition model 
(left) the deposit has an open structure with many holes. 
However, in contrast to fractal patterns, these holes do not 
exist at all length scales and the resulting patterns are not 
fractal. 














face and form a compact deposit.”° In the absence of sur- 
face diffusion the deposit has an open structure with many 
holes. However, in contrast to fractal patterns, these holes 
do not form at all length scales and the resulting patterns are 
not fractal. 


In the processes that we have investigated, we have found 
that in contrast to the structure itself which is not fractual, 
the surface can be described by fractal concepts. We have 
characterized the roughness of the surface by how widely it 
fluctuates from its mean value."® By investigating how this 
width varies with time in different model system we have de- 
veloped a description for the fractal properties of rough sur- 
faces. We find that surfaces must be characterized by two 
different fractal dimensions. One describing the variation of 
surface width with time (or average height) and another for 
describing correlations in space. Mandelbrot! has in- 
troduced the word self-affine to denote objects which need 
more than one fractal dimension for their description. Fig- 
ure 8, shows that scaling the time dependence of the surface 
roughness according to our scaling description leads to the 
collapse of all the data into a single “universal” curve, in- 
dicating that rough surfaces are self-affine fractals. Similar 
plots have also been obtained for a variety of other growth 
processes. These show that the fractal concepts are powerful 
tools for describing rough surfaces. 





Figure 8: 


Scaling the time dependence of the surface roughness 
according to our theory leads to the collapse of all the data 
into a single curve, indicating that rough surfaces are self- 
affine fractals. Similar plots have also been obtained for a 
variety of other growth processes. These show that the 
fractal concepts are powerful tools for describing rough 
surfaces. [From ref. 19] 
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Fractal Patterns in 
Diffusive Growth 


The formation of a wide variety of patterns is controlled 
by the strength of the gradient of a field at the interface be- 
tween the structure and the outside.** For example, the rate 
of growth in aerosol and colloidal aggregation depends on 
the concentration of diffusing particles, in solidification it 
depends on the temperature gradient, and in dielectric break- 
down it is proportional to the electric field. The common 
denominator among many growth processes is that the con- 
trolling fields, like temperature, concentration, pressure, and 
electrical potential, satisfy either the diffusion equation or 
the Laplace equation which is a special case of the diffusion 
equation. The differences in morphologies are often reflec- 
tions of variations in boundary conditions and external 
parameters. 
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The most characteristic feature of diffusive growth 
processes is instability. The reason goes back to the fact that 
this type of growth is governed by the gradient of a field. 
As in electrostatic, the gradient of the field is largest at sharp 
tips or points of high curvature. Thus, when a small piece 
of the front gets ahead of its neighbors, it gains an advantage 
because the gradient at that point will become larger and this 
tip will have a higher probability to grow. This instability 
forces the boundary to break-up into dendritic arms that grow 
faster and faster, enventually forming a tree-like fractal pat- 
tern. In some processes such as viscous fingering *’ and 
dendritic solidification,” * the surface tension tends to com- 
pete with this instability and prevent the patterns from be- 
coming completely random fractals. 





Figure 9: 


An example of a diffusion-limited aggregation [DLA] cluster. 
In DLA, particles are launched from far away and allowed 
to diffuse. Once a diffusing particle comes in contact with a 
growing aggregate it sticks to it and becomes part of it. 
This process of releasing random walkers continues until a 
large cluster is formed. The simulations indicate that DLA 
clusters of the type shown above are fractal with a 
dimensionality D~+1.7 in two dimensions and D~2.5 in three 
dimensions. 
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One of the consequences of the tip instability is that the 
mathematics of the problem is highly non-linear. In addition, 
it is subject to a boundary condition that is continuously mov- 
ing. Therefore, progress in solving the problem of diffusive 
growth processes from an analytic approach has been slow. 
On the other hand, in recent years numerical methods have 
been developed that are very effective in producing a variety 
of the patterns that are produced by diffusion-limited or Lap- 
lacian growth processes. The most important of thse methods 
is the random walk simulation technique used to study 
diffusion-limited aggregation [DLA]. In DLA, particles are 
launched from far away and allowed to diffuse. Once a diffus- 
ing particle comes in contact with a growing aggregate it sticks 
to it and becomes part of it. This process of releasing ran- 
dom walkers continues until a large cluster is formed. The 
simulations indicate that DLA clusters (see Figure 9) are frac- 
tal with a dimensionality D~1.7 in two dimensions and 
D~2.5 in three dimensions. Due to the importance of diffu- 
sive growth processes, the random walk approach has been 
used extensively in solving the diffusion equation in 
both fractal growth**-* and diffusion-limited pattern 
formation.*’-*° 


Viscous Fingering 


An example of the application of DLA model for study- 
ing diffusion-limited pattern formation is shown in Figure 
i0. The picture on the left shows an experiment and the one 
on the right is the simulation ®’ of a phenomenon known as 
viscous fingering. Viscous fingers are patterns that are formed 
when a less viscous fluid, such as water, pushes against a 
more viscous fluid, such as oil. The water, instead of push- 
ing the oil, tunnels through it in the shape of long thin fingers, 
leaving much of the oil behind. This practical problem, which 
is of considerable interest to the oil and energy industry, is 
an example of pattern formation in a Laplacian field. The 
reason is that the pressure at the interface between the two 
fluids satisfies the Laplace equation. Since the DLA model 
solves the diffusion equation, or the Laplace equation, it can 
be used to simulate the formation of viscous fingers. There 
is however an important difference between DLA and the two- 
fluid problem. In viscous fingering surface tension tends to 
smooth out the irregularities at the interface, whereas there 
is no surface tension in DLA. 

By incorporating surface tension and surface diffusion 
effects in DLA, we have effectively simulated the formation 
of viscous fingering patterns,”” as shown in Figure 10. We 
find in agreement with the experiments, that as the surface 
tension is decreased the patterns become more ramified. 
When the two fluids are miscible, the surface tension van- 
ishes and the resulting patterns are the same as ordinary DLA. 
This type of simulations and some of the recent experiments 
have provided new concepts and methods for investigating 
this complex and important problem. 








Figure 10: 


Simulations of viscous fingering using the diffusion-limited 
aggregations process. Viscous fingers are patterns that are 
formed when a less viscous fluid, such as water, pushes 
against a more viscous fluid, such as oil. The water, instead 
of pushing the oil, tunnels through it in the shape of long 
thin fingers, leaving much of the oil behind. As they grow, 
the fingers become narrower and narrower, until they re- 
duce to DLA type patterns (Fig. 9). [From ref. 27] 
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Laplacian Fractals 


An alternative to the random walk approach is to directly 
solve the Laplace equation on a lattice and determine the field 
everywhere. The growth probability is then determined from 
the gradient of the field.2? We have developed a model for 
studying growth in a Laplacian field which in analogy with 
the DLA model is based on an aggregation process.*° How- 
ever, in contrast to the DLA type models, more than one par- 


ticle is added to the growing cluster in each time step. The 
advantage of this approach is that it reduces the random fluc- 
tuations which are inherent in the random addition of single 
particles in DLA-like models. In our model, the field is cal- 
culated by solving the Laplace equation on a regular lattice 
subject to the constraint that on the cluster the field is zero 
and on a large circle around the cluster the field is constant. 
After the field is calculated, the gradients at all the surface 
sites are determined. These gradients are then normalized 
to the maximum gradient at the surface by dividing them by 
the value of the largest gradient on the interface. The nor- 
malized gradient at each site is compared with a parameter 
p and all perimeter sites with a gradient larger than p are 
filled. 

We first consider the case in which p is random. In the 
simulation,®® the value of the normalized gradient at a pe- 
rimeter site is compared with a random number and the site 
is occupied if the value of the gradient is larger than the ran- 
dom number. The process of picking a new random number 
for each surface site and testing and occupying the site is 
repeated until all the perimeter sites have been checked. A 
typical fractal pattern obtained from this type of simulation 
is shown in Figure 11. The most striking difference between 
Figures 9 and II is the highly anisotropic shape of the clus- 
ter in Figure ll. This four-fold anisotropy, which is due to 
the fact that the cluster was grown on a square lattice, may 
be used to simulate the effects of anisotropy due to crystal 
structure or surface tension which are present in many pat- 
tern forming processes. One of the important consequences 
of this anisotropy is that the fractal dimenstion of the cluster 
with anisotropy is smaller than the fractal dimension of an 
aggregate without anisotropy. Analytical *'** and numeri- 
cal*? studies indicate that in two dimensions the fractal 
dimension of an anisotropic cluster is about 1.5 rather than 
the value 1.7 for off-lattice DLA. This type of strong response 
to anisotropy is observed in many diffusion-limited 
phenomena, such as directional solidification and viscous fin- 
gering. It appears that anisotropy is an inherent property of 
many random fractal clusters.** 

If we let p be a fixed constant, then the model becomes 
a deterministic growth process,”* because in each growth 
step all perimeter sites having a normalized gradient larger 
than p are occupied. As a result, the patterns generated in 
this process are exact fractals, which we call Laplace Car- 
pets.*® Four examples are shown in Figure 12. In the limit 
p=0, all the perimeter sites are filled and the result is a dense 
polygon having the symmetry of the underlying lattice. As 
shown in Figure 12, for finite p, the patterns are regular frac- 
tals with a fractal dimension which decreases from 2 to 1 
as p is increased from 0 to 1. The most novel feature of the 
Laplace carpets is that the fractal dimension in a (determinis- 
tic) Laplace growth process can be tuned continuously. An- 
other important property of Laplace Carpets is their close 
similarity with cellular automaton. This is perhaps the first 
time that a process with long-range correlations has produced 
patterns similar to a process like a cellular automaton which 
is based on short-range, local rules. 
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Figure 11: 


A typical fractal pattern obtained by solving the Laplace 
equation directly on a square lattice. The most striking 
difference between this cluster and the DLA pattern (Fig. 9) 
is the highly anisotropic shape of this cluster. This 
anisotropy, which is due to the fact that the cluster was 
grown on a square lattice, may be used to simulate the 
effects of anisotropy due to crystal structure or surface 
tension which are present in many pattern forming 
processes (for example, see Fig. 13). [From ref. 30] 





Figure 12: 


Four examples of Laplace carpets. In the limit p=O, all the 
perimeter sites are filled and the result is a dense polygon 
having the symmetry of the underlying lattice. As can be 
seen from these patterns, for finite p, the patterns are 
regular fractals with a fractal dimension which decreases 
from 2 to 1 as p is increased from 0 to 1. 
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Dendritic Solidification: 
Snowflakes 


We have extended our model for growth in a Laplacian 
field to simulate dendritic solidification.** ** The growth of 
dendritic crystals is one of the most profound examples of 
pattern formation in a diffusion-limited growth process where 
the temperature is the relevant field. The formation of snow- 
flakes is perhaps the best known and puzzling example of 
these processes. Although subject to intensive efforts, previ- 
ous methods had not produced such regular and intricate den- 
dritic structures as are found in nature. Therefore, the de- 
velopment of the deterministic growth model ** which for 
the first time was capable of producing realistic dendritic pat- 
terns, including snowflakes, has been an important step in 
developing effective methods for understanding the forma- 
tion of snowflakes and dendritic crystals. 





One of the key ingredients needed for the description of 
dendritic growth is that, by heat conservation, the normal 
velocity at a point on the interface is proportional to the tem- 
perature gradient. In order to implement this boundary con- 
dition it is necessary to make p a time dependent variable. 
We have, therefore, discretized the time interval into a piece- 
wise function which allows us to vary p and the growth rate 
at the interface with time. Using this method we can produce 
practically all types of observed two-dimensional dendritic 
patterns, including snowflakes, by changing the surface ten- 
sion and the growth parameters. In Figure 13, we show only 
a few selected examples of snowflakes generated by varying 
the environmental conditions during the growth. We have also 
investigated the fractal properties of the snowflake patterns 
by studying the rate of growth of the tips. We find ** that in 
dendritic patterns the growth rate is similar to the growth rate 
in anisotropic, but random, DLA-type fractal patterns, which 
have a fractal dimension of about 1.5.***'** On the other 
hand, the dense snow crystals, such as hexagonal plates, are 
not fractal and have a dimension D~2. 





Figure 13: 


We show only a few selected examples of snowflakes 
generated by varying the environmental conditions during 
the growth. On the left, three exampies of snow crystal 
patterns are shown which were generated by the 
deterministic growth model. Three patterns of natural 
snowflakes are shown on the right. [From ref. 28] 
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Conclusions 


The wide variety of problems discussed in this brief re- 
view suggest that one of the most puzzling and challenging 
problems in physics is the inevitability of fractals. Fractal 
space and time have been found in most areas of science and 
technology. The advent of fractal geometry has provided 
scientists and engineers with fresh mathematical tools for in- 
vestigating complex phenomena that had been viewed intrac- 
table by generations of scientists. In recent years, important 
progress has been made by applications of fractal concepts 
to understand a host of problems from disordered materials 
and complex fluids to chaos and turbulence. There is every 
indication, however, that despite these remarkable advances 
much work remains to be carried out. 
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An overview of the theoretical and experimental aspects 
of the recently-discovered superconducting compcunds, which 
resulted from research supported by The Office of Naval Re- 
search is presented. This overview is divided into three sec- 
tions. In the first section a review of some of the theoretical 
and computational works is presented under the subsections 
entitled pairing mechanisms, electronic structure calculations 
and thermophysical properties. In the second section surface 
and interface chemistry issues related to the fabrication and 
use of high-temperature superconductors for high- 
performance applications are presented. Specific issues that 
are discussed include metallization and the formation of sta- 
ble ohmic contacts, and chemically-stable overlayers that are 
suitable for passivation, protection and encapsulation of su- 
perconducting material structures that can then be used un- 
der a wide range of environmental conditions. Lastly, issues 
are discussed that are related to each of the bulk high- 
temperature superconducting ceramic oxides which have 
received so much attention the past two years. These include 
TL,Ba,Ca,Cu,O,, with a critical temperature of 125 K, which 
is the current record. 


Theory 


Since the discovery of high-temperature superconduc- 
tivity in the CuO ceramics in 1986, there has been an explo- 
sion of theoretical works relating to these systems. It is not 
possible, in a short review of theoretical work like this, to 
examine the essential ideas of even a small fraction of this 
large body of literature. We select, quite arbitrarily, some of 
the work we are aware of and discuss the essential ideas of 
these. We apologize, in advance, to all the research workers 
whose work we have not reviewed. 

This theoretical overview is divided into three sections: 
(1) pairing mechanism, (2) electronic structure calculations 
and finally, (3) thermophysical properties. In the section on 
pairing mechanism we review the basic microscopic mechan- 
isms that might be responsible for producing superconduct- 
ing states. This is a fundamental question and the success 
of any theoretical model addressing it should be judged in 
relation to the experiments it explains and the phenomena 
it predicts. The second part, dealing with electronic struc- 
ture calculations, is quite important in that such calculations 
give detailed information about conduction processes, den- 
sity of states and anisotropies ink-space. Finally, in the sec- 
tion on thermophysical properties we review works relating 
to phonon dispersion and soft modes, dynamics of tetragonal- 
to-orthorhombic phase transitions and the temperature and 
concentration dependence of the stability of these phases. 
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Pairing Mechanism: 


In the conventional superconductors a pair of electrons 
of opposite spin and momentum form a bound state which 
leads to a coherent and highly correlated many-body super- 
conducting state. The attractive interaction between the elec- 
trons of the pair is mediated by lattice vibrations (phonons) 
and the electrons overcome their strong Coulomb repulsion 
by staying away from each other in time (retardation). The 
basic question in the new superconductors is—What is the 
pairing mechanism? For a review of some of the early pair- 
ing models proposed for the high-T, materials we refer the 
reader to the work by Rice.' 

One of the leading theories for the high-T, superconduc- 

tors is the resonating-valence-band (RVB) model proposed 
by Anderson? and his coworkers.* According to this model, 
the pairing mechanism is magnetic in origin and not of the 
conventional BCS type. The starting point of the RVB the- 
ory is a two-dimensional Hubbard model at half-filling with 
strong on-site Coulomb repulsion U and an attractive inter- 
site hopping energy t. Without oxygen doping (i.e., dopant 
concentration d = 0), the ground state of the above model 
is expected to be a long-range anti-ferromagnetic (AF) 
state,* but Anderson argues that frustration might favor a 
RVB state over an AF ground state. The basic idea of the 
RVB theory is that strong electron-electron correlations re- 
sult in a separation of charge degrees of freedom from spin 
degrees of freedom. At low doping (d # 0) and temperature, 
the quasi-particle excitations are believed to be “holons” (i.e. , 
charge-carrying spinless particles) and “spinons” (i.e. , spin-'2 
chargeless particles). 
Superconductivity is due to formation of a condensate con- 
sisting of holon pairs. Since Bose-Einstein condensation is 
not possible in a strictly two-dimensional system, the inter- 
planar couplings are important in giving rise to super- 
conductivity. 

The RVB theory has been amplified or modified in 
several ways,’ °’ the details of which we won't go into here. 
Rice and Wang ° have proposed a model in which supercon- 
ductivity is due to condensation of a pair of bosons which 
gives rise to quasi-particle excitation energies which are simi- 
lar to those of RVB theory, but different from BCS theory. 
Rice and Wang, however, favor a phonon interaction which 
mediates the attraction between the boson pair in question. 
Coffey and Cox ° have given a nice summary of the essen- 
tial points of the RVB theory in Section II of their paper. 

Another theory that has been proposed is the spin-bag 
mechanism of Schrieffer’s group.”° The starting point of their 
theory is the strong AF spin order in the neighborhood of 
the superconducting transition (T,) and two-dimensional 
spin-correlations over a length L (200 A) in the neighbor- 
hood of the Neel transition (Ty > >T,). Doping creates ex- 
tra holes in the system and these are assumed to be localized 
over a length / (< <L). The spin of the hole couples to the 
spin density of the AF background through an exchange in- 
teraction. Within a given domain a hole produces an effec- 
tive potential well or bag in its vicinity in which the hole is 
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self-consistently trapped. The pairing is due to an effective 
attractive interaction between two holes which overcomes the 
short range Coulomb repulsion. 

The importance of spins in producing a superconduct- 
ing state has also been emphasized in a model proposed by 
Emery." He has proposed that oxygen doping in the 
214-material creates holes at the oxygen sites and that there 
is a narrow band of oxygen holes which couple strongly to 
the local spin configuration of the Cu-sites. This strong cou- 
pling is responsible for an attractive interaction between the 
oxygen holes. 

We mention in passing that other mechanisms such as 
plasmons ” and excitons * have also been proposed as pos- 
sible candidates for being the condensate of the supercon- 
ducting state. 

Recently doubts have been expressed about the impor- 
tance of the magnetic origin of the pairing interaction. The 
bismuth superconductors * (BaKBiO,) which are closely 
related in structure to the earlier superconductors and which 
are free of copper provide counter examples to the mag- 
netic origin of superconductivity. Another counter exam- 
ple ® is the 123-material in which Cu is substituted 100% by 
Ag bringing T, down to 40 K. This has ied some people to 
propose that local structure might play a role in producing 
a superconducting state. A recent double-well model” of 
oxygen motion indicates trends in this direction; it has been 
shown that a strong T-dependent electron-phonon coupling 
parameter could produce T, = 100 K. 


Electronic Structure Calculations: 


The clectronic structure '*-*° of the undoped compounds is 
of immense interest in getting a clue to the origin of super- 
conductivity. These parent compounds are La,CuO, 
(denoted by 214), YBa,Cu,0, (denoted by 123) and 
Bi,Sr,CaCu,O, (denoted by 2212). The important question 
is: Why is the ground staie of these compounds an antifer- 
romagnetic (AF) insulator? The band structure calculations, 
so far performed within the local density approximation, give 
results which are in contrast to the experimental situation— 
they all produce a metallic ground state. The reason for this 
discrepancy is believed to be the strong electron correlations 
in the CuO planes which are not adequately taken into ac- 
count in a band picture. It would be interesting to have elec- 
tronic structure calculations which take these correlations into 
account. In this regard, spin-polarized band structure calcu- 
lations are expected to yield improved results over conven- 
tional band structure calculations. However, the results °° *° 
of spin-polarized band structure calculations do not seem to 
be definitive. 





We would like to conclude this section by stating some 
of the important conclusions that have emerged from the band 
structure calculations. In these calculations '*: ' the impor- 
tance of the two-dimensional nature of the CuO planes was 
emphasized. The Copper d(x? —y2) orbitals and the neigh- 
boring oxygen p(x,y) orbitals interact to produce bonding o 
and antibonding o* oribtals. Similar results *°-** were also 
obtained for 123- and 2212-compounds. Another impor- 
tant result is that the antibonding band was positioned 
closer to the Fermi energy E,;. In the new 2212- 
compounds, a pair of slightly filled Bi 6p bands provide 
additional carriers in the Bi-O planes.*! A remarkable fea- 
ture in these compounds is the charge separations between 
the two Bi-O planes.”* 


Thermophysical Properties: 


In this brief discussion of the thermophysical properties 
of the high-T, oxide superconductors we restrict ourselves to 
theoretical investigations of the lattice dynamics of these sys- 
tems. In particular, we review the work on phonon modes 
and oxygen vacancy ordering and their influence on the struc- 
tural transitions of the 214 and 123 superconductors. Knowl- 
edge of the phonon spectrum and its dependence on the 
oxygen distribution may prove to be of prime importance in 
elucidating the mechanism of high-T, superconductivity. 
Several scenerios of how electron-phonon interactions can 
lead to high transition temperatures have been proposed in- 
cluding oxygen motion in double wells,” interlayer cou- 
pling *° *’ and coupling to soft quasicyclic modes associated 
with underconstrained nearest-neighbor rearrangements.** 
At the 195th National Meeting of the American Chemical So- 
ciety held at Los Angeles, California, in September 1988, 
three papers were presented which discussed theoretically the 
maximum attainable transition temperature for the copper 
oxide superconductors. William A. Goddard III and cowor- 
kers, on the basis of quantum-chemical calculations involv- 
ing triplet magnon pairing, believe that the upper limit of the 
transition temperature may be as high as 230 to 250 K.*° 
Daniel C. Mattis *° predicts an upper limit of 180 K on the 
basis of arguments concerning the physical overlap of Wan- 
nier functions on neighboring atoms within each planar struc- 
ture. R. Bruce King *' also favors an upper limit of 180 K, 
basing his arguments on a graph theory analysis of the 
skeletal-bonding topology in the metal-oxide compounds. 


There have been a limited number of theoretical investi- 
gations on phonon frequencies and eigenvectors of the La, , 
(Ba, Sr),CuO,°*-** and YBa,Cu,O,, *° ** superconductors. 
Unscreened lattice dynamical models ** **, yielding only the 
bare phonon frequencies, gave fair agreement with ex- 
perimentally determined total phonon density of states, mean 
square atomic vibrational amplitudes and Debye temperatures. 
Weber ** has shown that the effect of screening, spectrum 
renormalization, due to conduction electrons gives rise to 
large Kohn anomalies near the Brillouin-zone boundary in- 
volving oxygen breathing modes. It was pointed out by 
Chaplot ** that an additional effect of renormalization is to 


hybridize the high-frequency modes, dominated by oxygen, 
with the medium-frequency, heavier-nuclei modes thus mak- 
ing necessary an effective-mass description for the analysis 
of the isotope effect. Cohen, et al.** employing the potential- 
induced breathing model were able to predict the tetragonal- 
to-orthorhombic phase transition in La,CuO, as arising from 
an instability, softening, in the B,, tilting mode of the 
I4/mmm tetragonal structure. The onset of this transition has 
been recently derived from symmetry principles *. 

As is now well established, YBa,Cu,O0,, undergoes a 
tetragonal-to-orthorhombic phase transition for 1 > 6 > 0 
which is a consequence of oxygen vacancy ordering in the 
Cu-O basal plane. Several approaches have been undertaken 
to model this transition and describe the phase diagram in 
(T,d)-space. Most treatments are of the 2-D lattice gas type, 
employing first and second nearest-neighbor interactions for 
oxygen and vacancies and solved by mean field tech- 
niques *°**, Another approach, which has been successful in 
predicting the microstructure of the multiphase region, uti- 
lizes the method of concentration waves **; here the oxide is 
treated as an interstitial compound of ordered oxygen atoms 
and vacancies on a simple lattice °°. The model of Mattis * 
deserves mention since it specifically accounts for the copper- 
oxygen bonds and thus enables predictions to be made con- 
cerning the dielectric response of the various phases. The pic- 
ture that emerges from these studies is that at high tempera- 
tures, T > 750 K, and/or stoichiometries 6 < 0.5 a tetragonal 
phase exists with random ordering of oxygen and vacancies 
in the Cu-O basal plane. For material with d = 0 the super- 
conducting orthorhombic phase is stable with oxygen (atoms) 
and vacancies forming alternate chains along the crystal b- 
direction. At low temperatures and intermediate stoichiome- 
tries phase separation occurs with micro regions of tetragonal 
phase, orthorhombic phase and a second cell-doubled or- 
thorhombic phase **. Phase transitions between the 
tetragonal and orthorhombic phases are of a second-order 
disorder-order type. 

As a final note we mention the work on establishing the 
temperature domain in which thermal fluctuations and criti- 
cal behavior dominate. Many of the theories described in this 
review are based on mean field techniques which become in- 
valid when fluctuations are important. Estimates of the criti- 
cal region, Ginzburg criterion, are in the range & © (T.-T)/T, 
~ 0.1-0.7 °*-54. However, it has been pointed out that the 
breakdown of near-field behavior is progressive **; the abil- 
ity of mean-field theory to predict non-universal quantities 
(prefactors, GL-parameters and T.) is lost within a region 
EP (EG) %, Brout criterion *°, which may cover most of the 
superconducting regime. 
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Surfaces and Interfaces 


The full integration of bulk and thin film high- 
temperature superconductors into existing and new technol- 
ogies of high commercial value appears limited by a number 
of surface and interface materials issues. Many of these is- 
sues can be stated in general terms because they are shared 
by each type of ceramic superconductor (2-1-4, 1-2-3, or 
2-1-2-2). Others are more specific to the material under 
study, e.g., the toxic character of the Tl 2-1-2-2 com- 
pound. As new materials with even higher critical tempera- 
tures are developed, analogous problems will be encountered. 
Hence, the knowledge base developed for one class of ce- 
ramic material may also apply to others. 

First, there are issues related to materials synthesis so 
that structures can be fabricated with predetermined shapes, 
sizes, and current carrying ability. These range from macro- 
scopic to microscopic. Second, there are challenges related 
to the fabrication of superconducting thin films on a variety 
of substrates, with Si being an obvious choice from the per- 
spective of microelectronic devices. Third, there are issues 
related to the formation of stable ohmic contacts, particularly 
for small samples and thin films. Fourth, there are problems 
related to the passivation, protection or encapsulation of small 
structures such as fibers or thin films, so that the supercon- 
ducting oxides can be used under a wide range of environ- 
ments. These and a great many other issues raise challeng- 
ing chemical questions. Spectacular progress has been made 
in the few months that the new superconductors have been 
in existence, and we can anticipate that rapid progress will 
be made in the near future. These efforts that address fun- 
damental issues will bring us one step closer to realizing 
breakthroughs in technologies of high commercial value. 

It might be thought that surface and interface issues 
should be separated from those involved in the synthesis of 
the superconductors themselves. This is certainly not the case, 
however, because many of the most exciting opportunities for 
these materials will be in high performance applications. In 
these applications, the size of the superconducting elements 
will be comparable to those of the other components. Scal- 
ing down or shrinking the size of a structure exacerbates prob- 
lems related to interfacial phenomena. Indeed, the challenges 
of forming, contacting, and protecting a superconducting strip 
1 ym wide and 0.1 ym thick point to the intermingling of a 
wide range of materials issues. 

To date, issues related to bulk synthesis and surface or 
interface characterization have also been intimately tied. This 
can be seen by noting that the starting point for surface or 
interface research is a well-characterized bulk material. Un- 
fortunately for the surface scientist, early samples did not 
spring completely characterized from the firing furnaces, to 
paraphrase the springing-forth of Minerva completely armed 
from the head of Jupiter. Instead, early samples were sin- 
tered, were 70-90% dense, and individual grains were clad 
with other phases °°’. Difficulties in characterizing these 
interfaces and identifying their intrinsic properties are 
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reflected in the literature of the last year for the 1-2-3 and 
2-1-4 materials, and the last few months for the 2- 1-2-2 
materials. Only recently have bulk samples of sufficient qual- 
ity been available so that fracturing could provide a clean sur- 
face °°. Today, it is relatively routine to find single crystals 
having dimensions of greater than 1 mm, and several small 
companies are preparing to sell single crystals as large a 1 
cm x | cm so that full characterization can be achieved. 


Early studies indicated that polycrystalline, sintered sam- 
ples of the 1-2-3 and 2-1-4 materials degraded rapidly 
upon exposure to a range of environments, including H,O, 
CO, CO,, O, and solvents.°*®° More recent work suggests 
that degradation is much slower and that some of the early 
problems were related to the intergranular phases (e.g. car- 
bonates or cuprates). Early work also showed that exposure 
of the 1-2-3’s to high-intensity ultraviolet and X-ray pho- 
ton beams produces substantial changes in the surfaces.® 
Again, recent work has indicated that most of the changes 
are related to the presence of second phases and were not 
intrinsic to the superconductors. At the same time, studies 
by Rosenberg and coworker ® pointed out some of the de- 
tails of photon stimulated desorption. Electron beams of high 
current or high energy also induce damage. Exposure of these 
materials to energetic ions, such as used in Ar sputtering, 
leads to surface modification, structural changes, and the loss 
of superconductivity.°° This damage points to the fragile 
character of these superconductors but also indicates that it 
may be possible to selectively alter portions of a thin film, 
for example writing nonsuperconducting lines on a supercon- 
ducting template. 

An issue that has arisen repeatedly has been the possi- 
bility of oxygen loss through the surface at room tempera- 
ture. To our knowledge, there is no clear evidence that oxy- 
gen is lost under static vacuum conditions. Instead, the 
exposure of freshly prepared surfaces to ultrahigh vacuum 
leads to the chemisorption of residual gases from the vacuum 
system. Recent work with single crystals has provided evi- 
dence for the rearrangement of surface atoms after cleaving 
in vacuum.®* This has been attributed to transgranular frac- 
ture and the exposure of energetically unfavorable surfaces. 
Indeed, such effects can be understood in terms of the highly 
anisotropic unit cell, but the restructuring does not neces- 
sarily result in oxygen loss. 

Attempts to form contacts to surfaces or to investigate 
the electrical properties of the high-temperature super- 
conductors have often been complicated by nonreproducibil- 
ity. This can be related to the chemical processes that occur 
at these surfaces. As detailed studies of representative inter- 
faces have shown, there is a strong tendency for reactive me- 
tals to leach oxygen from the superconductor to form new 
metal oxide bonding configurations.** * The result of this in- 
terfacial chemistry is a heterogeneous transition region be- 
tween the buried superconductor and the metal film. In par- 
ticular, a cross section through an interface based on the 
1-2-3, 2-1-4 or 2-1-2-2 superconductors would show 
the superconductor; a region where oxygen has been removed, 





where the structure is likely to be disrupted, and which is 
not superconducting; a region where the metal adatoms have 
formed electrically-resistive oxides; and the metal overlayer, 
possibly containing oxygen and dissociated superconductor 
atoms in solution and at the surface. Such an interface is 
shown schematically in Fig. 1. These interfaces are metasta- 
ble because thermal processing will enhance oxygen trans- 
port to the metal layer and will increase the amount of sub- 
strate disruption. Certainly, these interfaces would not form 
the ohmic contacts desired in device applications. Likewise, 
the cladding of superconducting filaments with copper sheaths 
does not seem propitious since there will be interfacial 
interactions and the formation of a nonsuperconducting layer. 
The scale of this disrupted region is at least 50 A for room 
temperature metal deposition and is likely to be much larger 
if the interface is processed at a higher temperature. 

Ohmic contacts to these superconductors can be formed 
with Ag and Au overlayers.** These metals show minimal 
tendency to form oxides, and their deposition does not seri- 
ously disrupt the superconductor. For contracts, then, these 
metals appear to be the materials of choice. Two caveats 
should be noted, however. First, processing may lead to Au 
clustering, as has been observed for Au layers on many sub- 
strates. Second, if the surface to be metallized has been ex- 
posed to the air, it will be coated with hydrocarbons, water 
vapor and the like, and its superconducting properties will 
most likely be compromised. Attempts using Ar-ion sputter- 
ing to remove these contaminants prior to metallization have 
been successful, but almost certainly at the expense of the 
structural integrity of the surface. It remains to be seen 
whether the amount of surface damage (and loss of super- 
conductivity) compromises the effective use of this approach 
to the fabrication of ohmic contacts for devices. 

The protection of these superconductors is also crucial 
if they are to be integrated with other technologies. Much 
less is known about such passivation issues, although attempts 
are presently under way to develop protective overlayers. Two 
types of materials that have been examined so far and appear 
promising are metal oxides and non-metallic insulators.® 
Metal oxides have been formed by the deposition of metal 
atoms in an activated oxygen atmosphere to form a metal 
oxide precursor that does not leach oxygen from the substrate. 
The oxide layer then serves as a diffusion barrier against 
oxygen loss or atomic intermixing. To date, efforts have been 
successful with the oxides of Al and Bi, and it is likely that 
other oxides will prove to be effective as passivation layers 
and/or diffusion barriers. Those investigated so far are 
insulators, but future work may lead to conducting oxides that 
can serve as both contacts and passivators. 

The second type of material that shows promise for pas- 
sivation and electrical isolation is CaF,.°% This large- 
bandgap ionic insulator has a high static dielectric constant, 
can be readily evaporated from thermal sources in molecu- 
lar form, and shows no tendency to modify the superconduc- 
tor. As such, it may be useful as a dielectric layer in advanced 
devices. 


We also note that both organic and inorganic polymers 
may be useful as encapsulants. In fact, fully protective coat- 
ings for high temperature superconductors may need to be 
developed to meet the specific demands of large- and small- 
scale applications. These protective coatings will likely in- 
volve a multilayer structure that would consist of thin- or 
thick-films of metals, low dielectric ceramics, and possibly 
organic and inorganic polymeric films. These multifunctional 
coatings would be effective as a diffusion barrier, would with- 
stand environmental stress and temperature cycling, and would 
maintain strong adhesion. Specific opportunities exist for the 
synthesis of precursor molecules, such as organometallics, 
for chemical vapor deposition (CVD), or sol/gel approaches 
for low temperature fabrication of these multilayer protec- 
tive coatings. 


Bulk Superconducting 
Ceramic Oxides 


Superconducting oxides have been known since 1964, but 
until recently the intermetallic compounds showed higher su- 
perconducting temperatures. In 1975 research scientists at E. 
I. DuPont de Nemours ® discovered superconductivity in 
the system BaPb,_,Bi-,O, with a T, of 13 K. The structure 
for the superconducting composition in this system is only 
slightly distorted from the ideal cubic perovskite structure. 
It is generally accepted that a disproportionation of the Bi(IV) 
occurs, namely, 2Bi(IV) (6s') > Bi(III) (6s?) + Bi(¥)(6s°) 
at approximately 30 percent Bi. Sleight found that the best 
superconductors were single phase prepared by quenching 
from a rather restricted single-phase region, and hence these 
phases are actually metastable materials. At equilibrium con- 
ditions, two phases with different values of x would exist; 
the phase with a lower value of x would be metallic and with 
a higher value of x would be a semiconductor. It is impor- 
tant to keep in mind that the actual assignment of formal va- 
lence states is a convenient way of electron accounting; the 
actual states include appreciable admixing of anion functions. 
The system BaPb,_,BixO, should be studied further since it 
contains compositions showing the highest T, for any super- 
conductor not containing a transition element. Recently, for 
example, Cava and Batlogg * have shown that Ba,K, BiO, 
gave a T. of almost 30 K, which is considerably higher than 
the 13 K reported for BaPb,, 75Bi,, O,. 

La,CuO, was reported by Longo and Raccah ® to show 
an orthorhombic distortion of the K,NiF, structure with 
a=5.363 A, b=5.409 A and c=13.17 A. It was also 
reported ® ” that La,CuO, has a variable concentration of 
anion vacancies which may be represented as La,CuOQ,_,. 
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Superconductivity has been reported for some preparations 
of La,CuO,. However, there appears to be some question as 
to the stoichiometry of these products since only a small por- 
tion of the material seems to exhibit superconductivity.” 

The extent of the anion vacancies has been recently re- 
examined,” and the magnitude of this deficiency is less than 
can be unambiguously ascertained by direct thermogravimet- 
ric analysis which has a limit of accuracy in x of 0.01 for the 
composition La,CuO,_,. However, significant shifts in the 
Neel temperatures confirm a small variation in anion vacancy 
concentrations. 

In the La,_,A,CuO, phases (A=Ca,Sr,Ba) the substitu- 
tion of the alkaline earth cation for the rare earth depresses 
the tetragonal-to-orthorhombic transition temperature. The 
transition disappears completely at x> 0.2, which is about 
the composition for which superconductivity is no longer ob- 
served. Compositions of La,_.A,CuO, can also be pre- 
pared”: * where A is Cd(II) or Pb(II). However, these 
phases are not superconducting, and it appears that the more 
basic divalent cations are necessary to allow Cu(II) to coexist 
with 0?-. The existence of Cu(I) and Cu(III) in 
La,,Sr CuO, is consistent with the ESR spectra, which 
shows the absence of square planar Cu(II), and the Pauli- 
paramagnetic behavior over the temperature range from 77 
to 300 K. Since the Pauli-paramagnetic behavior of 
La, .Sr,CuO, is consistent with delocalized electrons, this 
would also indicate a high probability for the existence of 
Cu(1)-Cu(IID) formed as a result of disproportionation of 
Cu(II). * Subramanian et al.”° have recently substituted both 
sodium and potassium into La,CuO,, giving rise to the com- 
positions La,_,A,CuO, (A = Na,K). However, only the so- 
dium substituted samples exhibited superconducting behavior. 

The compound Ba,YCu,O, shows a superconducting 
transition of +93 K and crystallizes as a defect perovskite. 
The unit cell of Ba, YCu,O, is orthorhombic (Pmmm) with 
a=3.8198(1) A, b=3.8849(1) A and c=11.6762(3) A. The struc- 
ture may be considered as an oxygen-deficient perovskite with 
tripled unit cells due to Ba-Y ordering along the c-axis. For 
Ba, YCu,0,, the oxygens occupy 7/9 of the anion sites. One 
third of the copper is in 4-fold coordination and 2/3 are five- 
fold coordinated. A reversible structural transformation oc- 
curs with changing oxygen stoichiometry going from or- 
thorhombic at x=7.0 to tetragonal at x=6.0. ” The value x = 
7.0 is achieved by annealing in oxygen at 400-500°C, and 
this composition shows the sharpest superconducting transi- 
tion. It was shown by Davison et al.” that these materials 
are readily attacked by water and carbon dioxide in air to 
produce carbonates. 

Recently Maeda et al.” reported that a superconduct- 
ing transition of 120 K was obtained in the Bi/Sr/Ca/Cu/O 
system. The structure was determined for the composition 
Bi,Sr,CaCu,O, by several laboratories.’? *! 


In most of the studies reported to date on the Bi/Sr/Cu/O 
system measurements were made on single crystals selected 
from multiphase products. The group at DuPont selected platy 
crystals having the composition Bi,Sr, ,Ca,Cu,O,,, 
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(0.9>x>0.4) which showed a T.~95 K. Crystals of 
Bi,Sr,,Ca,Cu,O,,, for x=0.5 gave orthorhombic cell cons- 
tants a =5.399 A, b=5.414 A, c=30.904 A.” These cell 
dimensions are consistent with the results of other investiga- 
tors * *', The structure consists of pairs of CuO, sheets in- 
terleaved by Ca(Sr), alternating with double bismuthoxide 
layers. Sunshine et al.*’ have indicated that the addition of 
Pb to this system raises the T, above 100 K. There are now 
three groups of superconducting oxides which contain the 
mixed Cu(II)-Cu(III) oxidation states, namely La,,A"CuO, 
where A"=Ba,Sr,Ca; RBa,Cu,O, where R is almost any lan- 
thanide; and Bi,Sr,,Ca,Cu,O,,,. 

Sheng and Herman have recently reported * on a high- 
temperature superconducting phase in the system 
TI/Ba/Ca/Cu/O. Two phases were identified by Hazen et 
al.*°, namely T1/Ba/CaCu,O, and T1,Ba,Ca,Cu,0. Sleight 
et al.”?** have also reported on the structure of 
T1,Ba,CaCu,O, as well as T1,Ba,CuO,. In addition, the su- 
perconductor T1,Ba,Ca,Cu,0,, has been prepared by the 
DuPont group ® and shows the highest T, of any known bulk 
superconductor, namely 125 K. 

A series of oxides with high T, values has now been 
studied for the type (A™O),A,"Ca, ,Cu,O,>,, where A(II]) is 
Bi or Tl, A(II) is Ba or Sr, and n is the number of Cu-O 
sheets stacked. To, date, n=3 is the maximum number of 
stacked Cu-O sheets examined consecutively. There appears 
to be a general trend whereby T, increases as n increases. 
Unfortunately, these phases involve rather complex ordering, 
crystals of the phases are grown in sealed gold tubes, and 
excess reactants are always present. The toxicity as well as 
volatility of thallium, coupled with problems in obtaining 
reasonable quantities of homogeneous single-phase material, 
presents a challenge to the synthetic chemist. It will also be 
interesting to see if these materials are truly more stable over 
time than the La,,A,CuO, or RBa,Cu,O, phases. 

When superconductivity in oxides was first discovered 
in 1964, it created little excitement because the T.’s were 
lower than chose found for intermetallic compounds. How- 
ever, today there exists a number of oxides having metal- 
oxygen o* bands °° 8” which give the highest T, values of all 
known superconductors. Despite optimistic reports from 
numerous laboratories, a number of fundamental chemical 
problems exist which make commercialization of these ox- 
ides particularly challenging. The following is a brief sum- 
mary of these difficulties. 

The best superconductors in the BaPb,,Bi,O, system 
are metastable materials prepared by quenching from elevated 
temperatures. These compounds become two-phase if they 
are annealed at lower temperatures. In fact, even in the single- 
phase high-T. region there will be a tendency for Pb and Bi 
segregation on a short-range basis. This is probably why it 
is difficult to obtain sharp transistions in this system, regard- 
less of the value of x. 

The compounds which crystallize with the K,NiF, 
structure (La,CuO,, etc.) also have serious chemical prob- 
lems associated with them. The compound La,CuO, can be 
non-stoichiometric in several different ways. Most of the at- 





tention has emphasized the oxygen stoichiometry, i.c., 
La,CuO,,. There is a tendency to form intergrowth phases * 
of the type La,CuO,nLaCuO,, i.e., perovskite LaCuO, in- 
tergrowth with K,NiF,-type La,CuQ,. 

For all of the copper-containing supereconductors where 
copper is found to be present with a formai valence of Cu(II), 
we are actually dealing with the reactivity of holes located 
in the n* band * *’. This results in instability of these cop- 
per compounds when exposed to moisture and/or air. It is 
highly unlikely that treating the surface of such bulk com- 
pounds with a protective coating will do much to prevent their 
destruction under most operating conditions. The high mo- 
bility of copper in ternary oxides, even at room temperature, 
will add to the reactivity of superconductor surfaces with air 
and moisture. Furthermore, non-superconducting phases can 
readily form with stoichiometries close to those of the su- 
perconducting phases. Hence, pure superconducting materials 
with the necessary properties will be difficult to prepare. 

Single crystal growth can be achieved from melts of many 
of the complex superconducting copper oxides. However, 
many are found to melt incongruently, as well as to form in- 
tergrowths of nonsuperconducting phases. Twinning and poor 
crystal growth are the rule rather than the exception. 
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ONR Leads in 
Optimization 


Neal Glassman 
Office of Naval Research 


In October 1986, the Office of Naval Research under- 
took an accelerated research initiative in projective transfor- 
mations for optimization. The initiative was motivated by the 
startling announcement by Narendra Karmarkar and Bell 
Laboratories in the fall of 1984 of a new algorithm for linear 
programming. 

The new initiative exemplifies the response to the 
challenging and urgent optimization problems found through- 
out the Navy, with its complex multi-echelon inventory sys- 
tems, the world’s largest warehouses, industral-type job shops, 
daily multi-criterion manpower assignment problems, en- 
gineering design appiications, and strategic and tactical plan- 
ning and allocation decisions of the greatest importance. Dur- 
ing the last three decades, through ONR, the Navy has been 
investing in and developing the mathematics of optimization 
as well as specialized algorithms, codes, and applications. 

ONR supports optimization principally through its 
mathematics division. Although several individual programs 
within the division contribute to the effort, numerical analy- 
sis and discrete mathematics in particular, the focus for op- 
timization is in the operations research program. 


24 Naval Research Reviews 


The new algorithm, by Karmarkar’s account, yielded run- 
ning times for some problems that were hundreds of times 
faster than those of the simplex method. In addition to this 
astounding computational performance, the algorithm was sig- 
nificant in two respects. First, it reached an optimal solution 
in a number of iterations that is a polynomial func- 
tion of problem size. Second, it operated in the interior of 
the linear programming (LP)-polytope. A polytope is the 
higher dimensional analog of a polygon, which represents 
the set of all feasible points. All of these attributes contributed 
to the widespread enthusiasm that greeted Karmarkar’s 
announcement. 

The simplex method of the solution of linear program- 
ming problems operates on the surface of the LP-polytope. 
The algorithm moves from vertex to vertex until it finds a 
vertex where the objective function is optimal. Roughly 
speaking, the complexity of the network of edges and ver- 
tices that defines the surface of the polytope grows exponen- 
tially with the size of the problem. For this reason, the sim- 
plex method is an exponential algorithm: in the worst case, 
the number of iterations required to reach an optimal vertex 
is an exponential function of the number of variables and con- 
straints. The insight that an algorithm operating in the in- 
terior of the polytope might avoid such complexity is not new. 
In fact, since the invention of the simplex method in the 
mid-1940’s, hundreds of papers implementing this idea have 
been written. None, however, was successful in yielding faster 
running times on problems of interest. 

Khachian’s algorithm, one such interior method, was an- 
nounced in 1978. The first algorithm for linear programming 
with a polynomial performance bound, it verified the cor- 
rectness of the insight that an interior method could avoid 
much of the complexity associated with the simplex method. 
Despite its theoretical superiority for worst-case problems, 
Khachian’s algorithm proved to be impractically slow on real- 
world problems. The difficulty with it and the other previ- 
ously developed interior methods lies in finding the direc- 
tion toward the optimum from an arbitrary interior point. LP- 
polytopes derived from large real-world problems tend to be 
long and narrow, like carrot sticks, not shaped like a cube. 
For this reason, it is very difficult to find an interior point 
near the center and even more difficult to find a direction 
that does not quickly hit the boundary. Karmarkar’s contri- 
bution was the discovery that the application of a suitable 
transformation, a projective transformation, to the n- 
dimensional space and to a specified interior point produces 
a new problem with desirable properties in a transformed 
space. In the transformed space, the specified interior point 
is near the center of gravity of the polytope. From this point, 
a direction can easily be determined that is likely to produce 
substantial improvement in the objective function. 

The accelerated research initiative that began in October 
1986 has begun to address the important theoretical issues 
raised by Karmarkar’s algorithm and other interior methods 
for optimization. It has also begun to overcome the obsta- 
cles to effective implementation. Over the next two years, 
it will apply these techniques to selected Navy problems. 





There are intimate connections between Karmarkar’s al- 
gorithm and several other areas of mathematics. Any itera- 
tive scheme can be viewed as a system of first-order differen- 
tial equations by considering infinitesimal steps in the 
direction of the next iterate. Karmarkar’s algorithm can thus 
be considered as a method of solving such a system. This 
analogy may lead to the construction of other algorithms 
based on higher order methods of solution. Karmarkar’s al- 
gorithm is a nonlinear algorithm in that the projective trans- 
formation mentioned earlier creates a nonlinear objective 
function. As such, it is one of the few nonlinear algorithms 
with a proven global worst-case convergence bound. Exten- 
sions that yield a polynomial algorithm for quadratic 
programming have already been constructed. The relation- 
ship of Karmarkar’s algorithm to other algorithms for non- 
linear programming are being pursued in this initiative. 

Since its conception 40 years ago, the simplex method 
has become an “engine” for optimization. A great many al- 
gorithmic structures have been built with this engine as the 
driving mechanism. To carry the analogy further, if the sim- 
plex engine is replaced by a Karmarkar engine, operating on 
different principles, these structures must be examined to en- 
sure their compatibility. One of the most important of them 
is the body of algorithms that perform optimization over net- 
works. These algorithms, together with specialized list- 
processing and modeling techniques, constitute one of the tri- 
umphs of operations research in the last decade. The trans- 
lation from matrices and pivots to the language of nodes, arcs, 
and trees has profoundly affected our ability to solve many 
large, complex problems that, on the surface, bear littic 
resemblance to networks. These translations are now being 
accomplished. 

The last five years have also seen great progress in an- 
other technology that uses the simplex engine, integer 
programming, in which the solutions must have integral 
values. Current methodology is based on the solution of se- 
quences of related linear programs. Depending on the tech- 
niques used, the elements in the sequence are related by ad- 
dition of extra constraints or by the placement of additional 
bounds on the variables. In either case, the efficiency of these 
methods is crucially dependent on the use of information 
generated in the solution of the preceding elements of the se- 
quence. Techniques are being developed for incorporating this 
information within interior methods. 


Returning to the engine analogy, the simplex method has 
been fine-tuned over the years to operate very efficiently on 
routine problems and still be responsive in stressful situa- 
tions. The performance of Karmarkar’s algorithm, which has 
not had the advantage of such extensive development, be- 
comes more impressive by comparison. Such comparisons 
may be misleading, however. Production LP codes are typi- 
cally run at the beginning of each planning period, with the 
structure of the model unchanged but with changes in the data. 
The solution for the current period is usually close to that 
for the preceding period. The simplex method is able to start 
from the last solution, thereby greatly reducing the number 
of iterations required. No such advanced starts are yet possi- 
ble with Karmarkar’s algorithmn. 

Sensitivity information concerning the optimal solution 
is often as important as the solution itself. The simplex 
method provides basic solutions, i.e., solutions that are at 
a vertex of the LP-polytope. This feature permits the compu- 
tation of marginal costs—changes in the objective function 
corresponding to small changes in resource constraints. In- 
terior methods, on the other hand, grow numerically unsta- 
ble as they approach the boundary and thus cannot be used 
to obtain such sensitivity information. We are developing new 
algorithms to produce this required sensitivity information 
while staying within the interior framework. 

The Navy shares large and complex planning problems 
with the commercial world. The research being supported 
under ONR’s accelerated research initiative is intended to 
make the solution of such optimization problems much 
cheaper. This, in turn, may result in an even greater number 
of problems from such diverse areas as pattern recognition, 
engineering design, inverse methods, and signal processing 
being formulated as optimizations, thereby extending mathe- 
matical programming into new areas. The research is also 
intended to reduce the time needed to reach optimal solu- 
tions. If these efforts are successful, the importance of linear 
programming for real-time applications may come to equal 
its current preeminence as a planning tool. 
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SILICON 
CARBIDE 
COMES 
OF 

AGE 
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Recent advances in the synthesis, 


growth, and processing of single crys- 
tal silicon carbide render it a 
promising material for advances in 
electronics and optical applications 
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Silicon and Carbon are among the most abundant 
materials on earth. While carbon is extensively used by na- 
ture in the molecular structure of most forms of life, both 
silicon and carbon are increasingly used in high technology 
applications. Silicon carbide (SiC) is the only intermediate 
phase in the Si-C equilibrium system at one atmosphere pres- 
sure. It does not occur in nature. The first reported synthesis 
was accidental —a result of attempts by Berzelius to make di- 
amond. Acheson began the commercial manufacture of SiC 
in the United States. He also gave it the trade name “carborun- 
dum” from the combination of the Latin words carbon and 
corundum. Silicon carbide is also a genezic name for several 
different materials that contain SiC and additional compo- 
nents and phases which have been consolidated into a solid 
body using different high temperature process routes. Col- 
lectively, these various materials possess a variety of micros- 
tructural relationships among the various phases which, in 
part, manifest control of the associated physical properties. 
In one or more of its polycrystalling forms, SiC is a prin- 
cipal candidate for high-temperature and/or high strength and 
and/or abrasion resistant applications, for instance, as heat 
exchangers in fossil fuel systems, critical parts for uncooled 
gas-turbine and adiabatic diesel engines and seals for pumps 
which move abrasive slurries. 

Silicon carbide is also a semiconductor. Active devices 
have been produced in its high purity (or controllably doped) 
single crystal form. However it is only now that some of these 
devices are becoming commercially available. Some of their 
potential applications are in high-temperature,-frequency, and 
-power electronic devices; blue-light emitting diodes; 
Schottky diodes; UV radiation detectors; high-temperature 
photocells; and heterojunction devices. It is this aspect of SiC 
which is the focus of this paper. 

Different crystal forms of the same chemical composi- 
tion are called polymorphs. Polymorphism commonly refers 
to a three-dimensional change affected by either a complete 
alteration of the crystal structure or a slight shift in bond an- 
gle. Polytypism is a special type of polymorphism which oc- 
curs in certain close-packed structures. In this phenomenon, 
two dimensions of the basic repeating unit cell remain con- 
stant for each crystal structure (polytype), while the third 
dimension is a variable integral multiple of a common unit 
perpendicular to the planes in the material having the highest 
density (closest packing) of atoms. Silicon carbide has long 
interested crystallographers as a primary example of poly- 
typism, for it has been determined to occur in more than 170 
different one-dimensional ordering sequences without appar- 
ent variation in stoichiometry (however, recent data suggests 
that the Si/C ratio decreases as a function of the degree of 
hexagonality in the crystal until an exact 1:1 ratio is obtained 
at the completely hexagonal (2H) polytype). 





The fundamental structural unit in all SiC polytypes is cova- 
lently bonded tetrahedron of four C atoms with a Si atom 
at the center. Each C is likewise surrounded by four Si atoms. 
The tetrahedra are linked through their corners, thus satisfy- 
ing the four-fold coordination at any strucutural point and 
occupying an array of positions analogous to those occupied 
by spheres in close-packed structures. The simplest polytype 
notation is due to Ransdell, who assigned the letters C, H 
and R to represent the cubic, hexagonal and rhombohedral 
structures, respectively (the only crystallographic structures 
found in SiC) and combined them with a numeral which 
represents the number of closest-packed layers in the repeat- 
ing sequency. The most common polytypes are 3C and 6H; 
4H, I5R and 2H have also been identified in crystallite form, 
but are much more rare. All other polytypes are combina- 
tions of these basic sequences. The lone cubic polytype (3C) 
is also referred to as beta-SiC; the remaining polytypes are 
known collectively as alpha-SiC. Although a considerable 
body of literature exists regarding the kinetics and ther- 
modynamics of polytype formation, growth, and stability, nei- 
ther the physiochemical mechanism that produces the peri- 
odic sequences nor the forces leading to structures with very 
large unit cells are understood for this material. This litera- 
ture does indicate that beta-SiC may be both the initial poly- 
type that forms at virtually all temperatures and a necessary 
precursor to the occurrence of other polytypes, as a result 
of transformation. This polytype is particularly interesting 
from the standpoint of electronics. 

A fundamental property of a semiconductor is the for- 
bidden region (bandgap) in which charge carrier energies are 
not allowed. Beta-SiC exhibits an electronic band gap of 2.2 
electron volts (eV) at 300°K. This is exactly twice that of the 
1.1 eV exhibited by the most common semiconductor (sili- 
con) and not closely approached by any other semiconduc- 
tor in common use. This large bandgap exists primarily as 
a result of the high strength of the chemical bond between 
silicon and carbon. This larger bandgap manifests itself in 


several ways. The most obvious one is that unlike common ' 


semiconductors, beta-SiC is transparent to light at wavelengths 
longer than 516 nanometers. Thus, while it absorbs blue and 
green light, it transmits red light very well. Between red and 
green in the spectrum lies yellow which is partially trans- 
mitted and partially reflected by beta-SiC. Figure | illustrates 
a film of beta silicon carbide lying over graph paper. The 
film is only 20 micrometes thick— about one half as thick 
as a human hair and was grown on and subsequently removed 
from a silicon substrate by the NASA Lewis Research Center. 





Figure 1 


Free-standing beta-Sic films courtesy NASA Lewis Research 
Center 














The larger bandgap exhibited by beta-SiC also translates 
to a much higher dielectric strength (electric field at which 
breakdown occurs). For SiC this number is 4 million volts/em 
or ten times that of the common gallium arsenide semicon- 
ductor and 13 times that of silicon. Since the output power 
which a semiconductor transistor is capable of delivering is 
proportional to the square of the voltage swing across that 
transistor, SiC is theoretically capable of delivering 100 times 
the output of a gallium arsenide device and 169 times that 
of a comparable silicon device. Another attribute of the higher 
bandgap is its thermal properties. As a semiconductor is 
heated, its electrons gain energy with the result that an in- 
creasingly larger number are transferred from the valence 
band to the conduction band. Since this number is inversely 
proportional to an exponential function of the bandgap, SiC 
transistors and diodes can be operated at a much higher tem- 
perature than can devices of common semiconductors. A 
striking example of this is shown in Figure 2. Here a beta- 
SiC transistor is shown under test at 650 degrees Celsius. 
At this temperature it glows red hot. This is fully 250 degrees 
hotter than any other transistor has ever been shown to oper- 
ate! The current versus voltage characteristics of the transis- 
tor are shown in Figure 3. Several aspects of SiC transister 
performance are notable. First, the curves have large horizon- 
tal sections which means that the output conductance is low—a 
very desirable property for an amplifier. Second—and most 
unique, is that the highest performance is obtained at tem- 
peratures in excess of 300°C. No other known semiconduc- 
tor increases its performance as a function of increasing tem- 
perature above its room temperature performance. Finally, 
performance at 650°C is nearly as good as room tempera- 
ture performance. This capability has rather profound im- 
plications. Aside from such down-to-earth applications as bet- 
ter control over (1) thermal appliances such as hair dryers, 
electric irons, and space heaters and (2) various functions 
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such as the air/fuel mixture in carburetor of automobiles, the 
space program could benefit substantially. Many satellites de- 
vote over 50% of their volume and weight to various forms 
of cooling equipment to cool their electronics to tempera- 
tures below that of boiling water. It can be shown that entire 
cooling systems in unmanned satellites can be eliminated if 
electronics can be found which will operate at temperatures 
in the vicinity of 350 degrees Celsius. Not only could the 
satellite be simplified with commensurate reliability improve- 
ments, but launching costs for the reduced payload weight 
could be substantially affected. 





Figure 2 
SiC transistor glowing red-hot at 650 C 








Related to the high temperature operation of electronics 
is the reliability of the semiconductors themselves. As tem- 
peratures increase, the impurities used to dope the semicon- 
ductors move or diffuse at proportionally higher rates. Since 
the bourdaries of these impurities are what define the tran- 
sistors and other devices made in the semiconductor, the 
properties of the semiconductor begin to change with impu- 
rity diffusion and the device can no longer perform its in- 
tended function. The extremely strong bond between silicon 
and carbon in the SiC crystalline lattice structure greatly in- 
hibits the diffusion of impurities. Figure 4 shows the diffu- 
sion of various impurities in SiC. The very impressive fea- 
ture is that at temperatures well above 1423 °C wherein the 
common silicon semiconductor melts, impurities in SiC are 
still rather stable! GaAs melts at an even lower temperature 
of 1220°C. Thus the typical diffusion dominated wear-out 
mechanism found in GaAs power transistors is not present 
in SiC at measurable rates. 
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The dielectric constant of a semiconductor helps to set 
the upper operating frequency of a device made from that 
material. In contrast to a dielectric constant of 11.8 for sili- 
con, 12.8 for gallium arsenide, and 14 for indium phosphide, 
beta-SiC exhibits a dielectric constant of only 9.72. This lower 
value significantly reduces capacitive parasitic losses as- 
sociated with the charging/discharging of transmission lines 
in circuits fabricated of SiC, and enables the reduction in size 
of the transistors needed to drive electrical signals down the 
microminiature transmission lines interconnecting the tran- 
sistors on chip and between chips. The decreased size fur- 
ther reduces power requirements and reduces power supply 
size and cost. 





Figure 3 
SiC transistor characteristics at 300 C and 650 C 
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The ability to dissipate heat is not only proportional to 
the temperature difference between the source and the sink 
(the efficacy of SiC to satellites), but is also proportional to 
the thermal conductivity of the material itself. Here SiC is 
also a winner. At room temperature its thermal conductivity 
is equivalent to that of copper (5 watts/em-°C). This com- 
pares to values of 1.5 and 0.5 for silicon and gallium arsenide, 
respectively. Thus while SiC can operate reliably at elevated 
temperatures (unlike common semiconductors), its high ther- 
mal conductivity also ensures that heat is adequately dissi- 
pated from the device. 








Figure 4 
Impurity Diffusion in SiC as a function of temperature 
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The velocity of the charge carriers (e.g., electrons and 
holes) in a semiconductor typically varies as a function of 
the electric field strength to which the carriers are subjected. 
In gallium arsenide and related III-V semiconductors this 
velocity can be very high for even comparatively low values 
of electric field. This engenders devices made of these 
materials to exhibit very high frequency response at low sig- 
nal levels. As signal levels increase and the required electric 
fields within the semiconductor increase proportionately, 
these materials lose their low field velocity advantages, as 
shown in Figure 5. Figure 6 provides a major part of the ex- 
planation for this phenomenon. Shown in red is the band di- 
agram for gallium arsenide. Gallium arsenide is a direct band- 
gap material in that the minimum of the conduction band 
occurs at the same point in “k” space (a manner of illustrat- 
ing momentum of charge carriers) as does the maximum of 
the valence band. As electrons in conduction band valley gain 
energy they are driven to the subsidiary conduction band val- 
ley where the parabolicity is less and their effective mass is 
greater. With the same energy, they must slow down in ac- 
cordance with the formula for kinetic energy. In contrast, SiC 
(like silicon) is an indirect bandgap material in that its con- 
duction band minimum does not lie in the same point in k 


space as does the valence band maxima. In SiC, the low 
parabolicity of the conduction band causes the electrons to 
be comparatively heavy and thus (for any given energy level) 
slower than electrons in gallium arsenide. Unlike electrons 
in gallium arsenide, however, electrons gaining energy in SiC 
can transfer to an upper subsidiary conduction band wherein 
the parabolicity is somewhat greater and the effective mass 
is somewhat lower thus causing them to gain velocity. In this 
respect, SiC is no different from silicon. As electrons gain 
further energy, however, they are scattered by their interac- 
tion with optical phonons. In both silicon and gallium ar- 
senide this scattering limits electrons to a velocity of 10 mil- 
lion cm/sec. In beta SiC, the lighter mass of carbon causes 
this limit to occur at a velocity 2.5 times greater. Thus while 
both silicon and gallium arsenide exhibit greater velocity at 
low electric field, SiC is clearly superior at high electric 
fields. 





Figure 5 
Velocity-field comparison of SiC with Silicon and GaAs 
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High charge carrier velocity at high electric fields is very 
important in two different types of electron device applica- 
tions. The first is in power devices operating in the micro- 
wave and millimeter wave spectrum. To obtain high power 
at microwave frequencies, the high carrier velocity must be 
present at high electric fields. The very high electric field 
breakdown (dielectric strength) of SiC enables large voltage 
changes simultaneously with high carrier velocity. Figure 5 
illustrates the order of magnitude superiority of SiC in this 
respect. To be able to operate at frequencies in the millimeter 
wave spectrum, the velocity of the electrons must be high. 
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Figure 6 
Band Diagrams for GaAs and SiC 
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It is seen in Figure 5 that electrons in silicon can never reach 
the velocity-field product required for microwave power am- 
plifiers (2 x 10’ cm/sec and 5 X 10‘ volts/cm). While gal- 
lium arsenide meets the velocity requirements for low elec- 
tric field strengths, it falls short at the higher fields needed 
for high power amplifiers or oscillators. Over two decades 
ago and prior to any serious attempts to synthesize semicon- 
ducting grade SiC, A. Johnson published a method for select- 
ing the proper semiconducting material for high performance 
electron devices. Table I is derived from those criteria and 
compares various candidate semiconductor materials. The 
potential of SiC is striking! 
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The second type of electron device wherein high charge 
carrier velocity at high electric field strength is clearly an 
attribute is that of high density integrated circuits. As the den- 
sity of logic functions per unit area increases, individual tran- 
sistor size must decrease as indeed it has done over the past 
four decades. Decreased transistor size also provides a shorter 
path length (channel) for the electrons to transit between the 
transistor source and drain as shown in Figure 7a. Owing 
primarily to noise margin limitations, logic swing levels (volt- 
age levels between a logic 0 and logic 1) must be maintained 
at several volts. (Typical logic levels today vary between 0 
and 5 volts.) With the applied voltage fixed and the transis- 
tor channel length decreasing, the internal electric fields con- 
tinue to increase. State-of-the-art transistors exhibit one 
micrometer channel lengths. At this value it is seen from Fig- 
ure 5 that both silicon and gallium arsenide exhibit saturated 
velocities of but 10 million cm/sec whereas SiC still exhibits 
a positive mobility (velocity versus field slope) and is com- 
fortably far from its breakdown limits. In 1972, R.W. Keyes 
developed a figure-of-merit for the selection of a semicon- 
ductor material for use in dense integrated circuits. Factors 
used by Keyes include dielectric constant, thermal conduc- 
tivity and saturated velocity—in all of which SiC excels. Ta- 
ble II has been compiled based on the criteria of Keyes. From 
it, there are strong reasons to use SiC in integrated circuit 
development. 

The low field mobility (velocity versus field slope) of 
silicon carbide is lower than that of common semiconduc- 
tors as shown in Figure 5. While that could be interpreted 
as a serious disadvantage for devices manufactured only 5 
years ago, present generation micrometer-dimensioned elec- 
tron devices are much more dependent on saturated velocity 
than on low field mobility. Even newer device configurations 
may provide a decided advantage to SiC. One such configu- 
ration is shown in Figure 7b. This device is known as a dual 
gate field effect transistor (FET) wherein electrons transit 
from the source to the drain and are modulated by the gates. 





Table I. 
Johnson Figure of Merit 





E(BD) V(sat) 


volts/cm cm/sec 


MATERIAL 
Silicon 
GaAs 

InP 

SiC (alpha) 
SiC (beta) 


E(BD)*V(sat) 


Pwr*Freq? *imped 
Tl 
volts/sec Watt-Ohms 
Sec? 
x10" x1028 

95 9.0 
25.0 62.5 
38.0 144.4 

250.0 6250.0 
320.0 10240.0 





Table Il. 
Keyes Figure of Merit 





Material 


Units 


Silicon 
GaAs 

InP 

SiC (beta) 


K 


dimensionless 


11.8 
12.8 
14.0 

9.7 





In theory, the dual gate device provides much higher gain 
and much better output conductance than does a single gate 
device. When fabricated in gallium arsenide, electrons un- 
der the second gate and beyond are in their heavy mass, lower 
velocity state as their energy levels are high. As such, the 
second gate cannot be used to modulate the device at micro- 
wave frequencies. In SiC, however, electrons transiting un- 
der the second gate are at their very fastest velocity. In fact, 
the AVERAGE velocity of the electrons in the SiC device is 
expected to be substantially higher than in the gallium ar- 
senide device. In this respect, the dual gate FET effectively 
circumvents the lower low field mobility of SiC. If the first 
portion of this dual gate SiC FET is used solely as an elec- 
tron injector for the second part of the device, and the over- 
all device is submicrometer so as to further reduce the prob- 
ability of electron scattering, then not only will the average 
velocity of electrons be much greater than that of electrons 
in most common semiconductors, but the electrons will all 
be traversing at velocities nearly equal to each other. If such 
can be verified, SiC (contrary to intuition) could also exhibit 
a very low noise background and thus compete with the now- 
dominant gallium arsenide devices in this application. 


Other reasons for choosing SiC are related primarily to 
immunity to radiation damage and to burn out from elec- 
tromagnetic pulses (EMP) such as those originated by light- 
ning. While no measurements have yet been conducted to 
verify the immunity of SiC circuits to EMP damage, the in- 
trinsic advantages of higher dielectric strength and higher 
thermal conductivity will certainly be beneficial. Recent co- 
balt 60 gamma ray measurements have shown SiC field ef- 
fect transistors to be unaffected to total doses exceeding 10’ 
RADS. 
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Figure 7 


(a) Single gate Field Effect Transistor 
(b) Dual gate Field Effect Transistor 
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Optical Properties. Microparticle abrasion is a serious 
problem in many optical lenses, mirrors, and windows. The 
hardness of SiC makes it potentially useful in airborne, 
spaceborne, and related applications. The optical reflectivity 
of SiC is shown in Figure 8. Although it has strong reflec- 
tivity at 12 micrometers, through the optical spectrum of 0.6 
to 10 micrometers, it is quite acceptable for use as a window. 
Undoped SiC in its alpha form is transparent and colorless, 
but aluminum doping turns it blue and boron doping colors 
it brown. Beta-SiC (as was earlier mentioned) is yellow 
because of its bandgap, but the addition of nitrogen impurities 
tends to color it green. Hexagonal 6H material that is nitrogen 
doped is green; whereas, undoped 6H material is colorless. 
The index of refraction of the cubic beta material is 2.48 while 
that of the hexagonal material tends toward 2.7. 

Nucleation, Synthesis, and Growth. Many of the above- 
stated properties and applications of SiC have long been 
known, but the required SiC materials simply could not be 
produced in the required perfection, quantities, and yields. 
During 1987/8 most of these limitations have been overcome. 
The most significant of these was the synthesis and repro- 
ducible growth of single crystal SiC boules of the 6H hex- 
agonal type via sublimation growth techniques. Another im- 
portant factor in the growth of usable SiC boules was the 
improvement in purity of the starting materials. A porous cyl- 
inder of graphite is coaxially located inside an outer, non- 
porous graphite clyinder. Between the two is placed a mix- 
ture of high purity powdered silicon and carbon. The seed 
crystal is placed inside the porous cylinder and the entire ap- 
paratus is brought to a temperature of 2250°C with a 
temperature gradient in the vicinity of the seed which is main- 
tained as the coolest part of the apparatus. Monocrystalline 
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Figure 8 
Optical Reflectivity of SiC 
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boules of 6H SiC grown from this apparatus have exhibited 
crystalline defect densities as low as 10*/cm*—similar to 
that of gallium arsenide crystals used today. Figure 9 depicts 
one such boule. Although current boule diameters are no 
larger than 2 cm in diameter, 2 inch diameter boules are under 
active development. Even then, the 2 cm diameter boules pro- 
vide material adequate in size to serve as windows and suitable 
for use in electron device development. 

While the process has not yet been perfected to grow SiC 
boules of cubic crystalline structure, research directed toward 
that pursuit is currently underway. Cubic SiC films have been 
epitaxially grown on the hexagonal substrates. When the 6H 
boules are cut into wafers such that the basal plane (c-axis) 
is the wafer surface, cubic SiC is grown by chemical vapor 
deposition methods flowing silane and ethylene gases car- 
ried in H, over the wafer at temperatures in the vicinity of 
1400°C. In this juxtaposition orientation the crystallographic 
axes of the growing cubic and the underlying hexagonal 
crystals perfectly align as shown in Figure 10. Here the in- 
terface between the hexagonal substrate and the cubic epitaxial 
layer is viewed using high energy, high resolution transmis- 
sion electron microscopy. The layered hexagonal underlying 
crystalline structure is easily discernable from the overlying 
epitaxial cubic film. The atoms of each crystal are seen to 
be perfectly aligned at the interface. Should the surface of 
the underlying 6H hexagonal wafer be misoriented more than 
2 degrees off the basal plane, only hexagonal epitaxial films 
of 6H—SiC can be grown. 








Summary. While SiC is clearly superior to the conven- 
tional semiconductors in many critical respects and it has at 
6H SiC crystal boule last (in small size) been produced with perfection as good 
as what is now commonly used in the gallium arsenide in- 
dustry, much remains to be done. The investment in SiC to 
date is less than 1% of the total investment in gallium arsenide 
which, in turn, has only experienced a few percent of the 
investment experienced in silicon. Although processing of a 
material with the hardness of SiC is difficult, it is not with- 
out advantages in that it is not easily damaged during 
manufacture. In this respect industrial yield may ultimately 
exceed that of gallium arsenide. If past continues to be pro- 
logue we can expect that the established semiconductor in- 
dustry will not espouse SiC. Just as the established ger- 
manium industry left newcomers to exploit silicon, and the 
silicon industry left newcomers to exploit gallium arsenide, 
the SiC industry is expected to grow from newly-formed in- 
dustries championing its cause. Even then, we can expect to 
see the first SiC electronic products in everyday use early 
in the next decade. 


Figure 9 















































Figure 10 Biographies 


Cubic SiC on hexagonal SiC Interface Max Yoder is an electronics engineer with the 
Solid State Electronics Branch of the Office of Naval 
Research (ONR) where he supervises programs in 
semiconductor materials, monolithic microwave inte- 
grated circuits, gigabit logic, new electronic device 
structures, and electromagnetic system applications. 
He has published and presented many papers in these 
fields of solid state materials. He is a pioneer in the 
development of research programs for monolithic 
microwave devices in the U.S. Currently he manages at 
ONR the research program in diamond semiconductor 
technology. 





Dr. Robert Davis is Professor of Materials Sci- 
ence and Engineering at North Carolina State Univer- 
sity. His research and papers have been in the fields of 
thin films for microelectronic devices, self-diffusion 
of 14C and 39Si. and high temperature creep in silicon 
carbide. He is a Fellow of the American Ceramic 
Society and recipient of the ALCOA Distinguished 
Research Award. 
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Research Notes 


Navy Honors 

Dr. Fred Saalfeld 
for 

Science Research 
Management 


The Honorable Thomas F. Faught, Jr., (left) Assis- 
tant Secretary of the Navy for Research, Engineering 
and Systems, presented the Conrad Award for 1988 to 
Dr. Fred E. Saalfeld, Director of the Office of Naval 
Research, who is accompanied by Mrs. Saalfeld. 


Dr. Fred E. Saalfeld, Director of the Office of 
Naval Research (ONR), is this year’s recipient of the 
U.S. Navy’s Robert Dexter Conrad Award, the highest 
Navy honor for scientific achievement. 

The award, in the form of a gold medal and a 
certificate, was presented to Dr. Saalfeld by Thomas F. 
Faught, Jr., Assistant Secretary of the Navy for Re- 
search, Engineering and Systems, during a ceremony in 
the Office of the Chief of Naval Research. 

The Conrad Award is presented annually to an 
individual who has made outstanding contributions in 
the field of research and development for the De- 
partment of the Navy. It commemorates Captain 
Robert Dexter Conrad who was one of the primary 
architects of the Office of Naval Research. the R&D 
organization responsible for the basic research pro- 
grams of the Navy. 

Dr. Saalfeld has been working for the Navy since 
1962 when he joined the Naval Research Laboratory 
(NRL), a field activity of the Office of Chief of Naval 
Research. At NRL, he worked as a bench chemist, 
head of the Mass Spectrometry Section, then as head 
of the Physical Chemistry Branch. In 1976, he became 
Superintendent of the Chemistry Division. During his 
tenure as Superintendent, he also served as Chief 
Scientist of ONR’s Branch Office in London. 

Dr. Saalfeld was appointed Director of the ONR 
Contract Research Department in 1982 and Associate 
Director of ONR in 1985. In 1987, he became Director 
of the Office of Naval Research . 

Mr. Faught praised Dr. Saalfeld for his excep- 
tional contributions and outstanding career as a sci- 
entist and manager. In presenting the award. Mr. 
Faught said of Dr. Saalfeld’s work: 

“... As a bench scientist and science manager in 
the Naval Research Laboratory’s chemistry program. 
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Dr. Saalfeld developed the Central Atmosphere 
Monitoring System (CAMS), now in service on the 
Navy’s nuclear submarines, and was instrumental in 
developing new secondary ion mass spectroscopy 
techniques that are widely used today in molecular 
biology and medical research. As Direcfor of the Of- 
fice of Naval Research, Dr. Saalfeld has aggressively 
stressed and strengthened the scientific excellence of 
the Navy’s Basic Research Program. His innovative 
management initiatives, exceptional scientific leader- 
ship, and outstanding management accomplishments 
have confirmed the position of the Office of Naval 
Research as a leader among federal agencies. His dis- 
tinguished scientific achievements, exemplary man- 
agement of research programs, and commitments to 
scientific excellence have directly contributed to the 
development of major new operating capabilities for 
the United States Navy.” 


Dr. Saalfeld is the 33rd recipient of the Conrad 
Award. Other awards and honors he has received 
include the Presidential rank of Meritorious Execu- 
tive, the Navy Superior and Meritorious Civilian Ser- 
vice Awards, and the Navy Superior Merit Award. 

Dr. Saalfeld is Founder and First President of the 
Mass Spectrometry Discussion Group of Greater 
Washington, Past National Councilor of the Ameri- 
can Chemical Society, Past Secretary of the American 
Society for Mass Spectrometry, and member of 14 
national and international professional societies. 

He has authored over 450 scientific papers. re- 
ports and presentations and has patents for three sci- 
entific inventions. 

A native of St. Louis, MO, Dr. Saalfeld obtained 
his Ph.D. from Iowa State University, and his B.S. 
from Southeast Missouri State University from which 
he received the Alumni Merit award in 1988. 








Research Notes 


Superconductors for the Navy 


Superconductivity, the transport of electrical current 
without energy loss due to resistance, and the devices that 
use this phenomenon are now leading topics in scientific and 
technological circles, including the Navy. The Navy never 
ignored the possibilities of superconductivity: quite to the con- 
trary. For example, the Office of Naval Research (ONR), 
since its founding in 1946, has supported superconductivity 
research. 

That year, ONR sponsored two conferences for scien- 
tists it supported at universities and in industry who were 
working in the field of superconductivity. From this meager 
beginning, ONR became a focal point for much of the low 
temperature physics and superconductivity research in the 
United States. Today, it manages some $22 million in Navy 
and Defense funds in this area. 

The interest today stems from the discovery, in 1986, of 
high-temperature superconductors. This important 
breakthrough occurred when scientists at IBM Zurich 
reported a material which became superconducting at 
temperatures as high as 30K(-243C). Soon after scientists 
at the University of Alabama and Houston discovered a ma- 
terial which superconducts at 90K(-183C). 

These new materials give hope to such forecasts as: trains 
levitated by magnetic forces: powerful, small and cheap com- 
puters; techniques for reading brain cell signals; electric cars; 
and electric motors which are half the size, weight, and cost 
of motors today. 

The phenomena of superconductivity was first recognized 
in 1911 when the Dutch physicist Kamerlingh Onnes noted 
that as he cooled mercury to the temperature of liquid helium 
(4K) the electrical resistance of the mercury dropped sud- 
denly from a finite value to zero. He concluded that mer- 
cury had passed on to a new state, which because of its 
extraordinary electrical property, he called the superconduc- 
ting state. Since then scientists have discovered that upon cool- 
ing, the zero electrical resistance state can be produced in 
about 1,000 elements, alloys and compounds, including such 
common metals as lead and tin, as well as a few exotic organic 
compounds. The temperature at which material transforms 
from its normal resistance state to a state of perfect electrical 
conductivity is called its transition temperature, T.. 

The transition temperature of most superconductors is 
very low, and because of this, the practical applications of 
superconductivity have been extremely limited. Over the 
years, the T. of newly discovered superconductors has risen 
very slowly, until recently the highest known T. was only 
23K in the compound Nb,Ge. Then in 1986, came the news 
of high-temperature superconductors producing a quantum 
jump in T, that will bring revolutionary results to the field, 
according to Dr. Donald Gubser who heads the supercon- 
ductivity program at the Naval Research Laboratory (NRL) 
in Washington, DC. 

For the Navy, ONR’s long term commitment to low 
temperature physics has already helped produce ultrasensitive, 
low noise, magnetic and electromagnetic detectors, electric 


drive systems currently under evaluation for ship propulsion, 
and ultra high speed electronic logic and memory circuits, 
all based on superconductivity. Civilian applications include 
magnetocardiography, magnetoencephalography, and some 
600 large supermagnets which are essential components of 
magnetic resonance imaging systems which are revolutioniz- 
ing noninvasive medical diagnosis. Other applications include 
the more than 1,000 supermagnets that guide protons around 
the four mile circumference of the Fermi Laboratory’s parti- 
cle accelerator and enormous supermagnets for the confine- 
ment of plasmas to achieve controlled thermonuclear fusion. 
ONR now supports research programs at NRL for 
developing the new high T, superconducting materials dis- 
covered in 1986 and a variety of processing techniques to 
further improve their properties. These new materials are 
made by mixing oxides and carbonates of yttrium, barivm, 
and copper, pressing into shape, and sintering at 1000C in 
a platinum crucible to form bulk samples. The samples are 
then tested for superconductivity by cooling in a special dewar 
where the electrical resistance is measured as a function of 
temperature. These samples are used to test a variety of pro- 
perties and to understand the materials from a microscope 
point of view. According to Dr. Gubser, the arrangement of 
atoms in the materials as well as their high transition 
temperature suggests that new mechanisms may produce a 
host of novel physical properties in these materials. He claims 
that “scientists must reexamine many of their long held con- 
cepts about the mechanisms producing superconductivity as 
well as the materials where it is most likely to occur.” 
From a technological viewpoint, the new materials will 
permit applications of superconductivity with much simplified 
refrigeration systems. Many superconductivity applications 
will now become technologically and economically viable 
options for a variety of uses. And in those cases where super- 
conductivity is the only available option for an emerging 
technology, these new materials will speed the development 
due to the significant easing of the refrigeration requirements. 
These superconductors, with critical transition 
temperatures above 77K (typical of a liquid nitrogen cryogen 
or a space environment), have caused the Navy RDT&E com- 
munity to reconsider use of superconducting components and 
devices. Potential areas for applications of these would 
include, antisubmarine warfare; antisurface warfare; antiair- 
craft warfare; and command, control, communications 
intelligence. Here they may: (a) perform the same function 
as previously known superconducting devices but at a more 
easily achievable temperature, (b) perform functions of 
devices which used older, low T, materials as above but 
better, or (c) perform entirely new functions, thus opening 
a broad array of opportunities for naval warfare systems. 
The U.S. Navy has capitalized on its many years of re- 
search in superconductivity and materials in rapidly respond- 
ing to these opportunities. Its current program is directed at 
both a basic understanding of these new materials and their 
application. 
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Tech Base Notes 


Described below is a synopsis of recent technol- 
ogy accomplishments sponsored by the Office of 
Naval Technology, which is the Navy research and 
development (R&D) organization responsible for 
conducting the exploratory development (6.2 catego- 
ry) for the Navy. Exploratory development programs 
are conducted in industry, academia and in the 
Navy’s own R&D centers. 


Corrosion Protection for Aircraft Components 


Naval aircraft structures that are subjected to high 
temperatures experience a high degree of corrosion 
and require frequent replacements of expensive com- 
ponents. Under Office of Naval Technology (ONT) 
sponsorship, the Naval Air Development Center has 
developed an organic coating that is thermally stable 
up to 700° F and provides corrosion protection far 
superior than materials currently used in military air- 
craft. This coating air dries at room temperature and is 
compatible with both steel components which are sub- 
jected to extreme heat. The new material consists of a 
protection by both the formation of a barrier layer and 
a chemical reaction on the substrate. The material is 
currently being tested by the Naval Air Depot A-4 and 
SH-60 aircraft, In addition, samples are being evaluat- 
ed by General Dynamics and McDonnell Aircraft 
Company for use on various other Navy aircraft. 
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Reducing Shipboard RF Radiation Hazards 


In the confines of a shipboard environment, RF 
transmissions can cause severe health-related effects to 
personnel as well as hazards of electromagnetic radia- 
tion to ordnance (HERO). Under Office of Naval 
Technology (ONT) sponsorship, the Naval Ocean Sys- 
tems Center has developed a method to calculate near 
field RF radiation for shipboard antennae which will 
be useful in reducing hazards. Digitized input of 
geometries from blue prints, isometric display of ship 
system models and graphics-mode editing of isometric 
images used in conjunction with this calculation 
method will expedite the assessment and redefining of 
operation practices and the redesign of modification of 
antenna site plans. These reconsiderations of restricted 
area designations and antennae locations will prevent 
certain health hazards and minimize HERO unsafe 
conditions. The American National Standards Insti- 
tute standard recently adopted by the Navy is much 
more restrictive than the previous standard, especially 
for the HF (2 to 30 MHz) band. This efficient, accurate 
prediction capability and anticipated long-term medi- 
cal cost savings make the method cost effective; the 
Naval Sea Systems command is using it to assure com- 
pliance with the new standard. 
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